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Abstract. We propose microwave-controlled rotations for qubits realized as Majorana bound 
states. To this end we study an inhomogeneous Kitaev chain in a microwave cavity. The chain 
consists of two topologically nontrivial regions separated by a topologically trivial, gapped 
region. The Majorana bound states at the interfaces between the left (right) regions and the 
central region are coupled, and their energies are split by virtual cotunneling processes. The 
amplitude for these cotunneling processes decreases exponentially in the number of sites of 
the gapped region, and the decay length diverges as the gap of the topologically trivial region 
closes. We demonstrate that microwave radiation can exponentially enhance the coupling 
between the Majorana bound states, both for classical and quantized electric fields. By 
solving the appropriate Liouville equation numerically we show that microwaves can drive 
Rabi oscillations in the Majorana sector. Our model emerges as an effective description for a 
topological semiconductor nanowire in a microwave cavity. Thus, our proposal provides an 
experimentally feasible way to obtain full single-qubit control necessary for universal quantum 
computation with Majorana qubits. 
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1. Introduction 

Majorana bound states (MBS) [1] are currently a strong focus of research in the condensed- 
matter community [2, 3, 4]. Semiconductor nanowires with strong Rashba spin-orbit coupling 
have emerged as a promising platform to host MBS. Following theoretical proposals [5, 6, 7] 
first experimental signatures of MBS have recently been reported [8, 9, 10, 11]. 

On the one hand, MBS are fascinating quantum systems in their own right. Despite 
the fact that Majorana fermions were predicted theoretically 70 years ago [12], it remains 
unclear to this date whether they exist as elementary particles. This means MBS may be 
their closest living relative. On the other hand, MBS could become useful for quantum 
computing because braiding MBS makes it possible to implement topologically-protected 
qubit operations which are immune against certain types of noise [13]. However, these 
braiding operations do not form a universal set of gates needed for quantum computation, 
so they have to be supplemented by other gates which are not topologically protected. 
The circuit quantum-electrodynamics (circuit-QED) architecture [14] offers a controlled 
and well-developed toolbox in the microwave domain and is thus an ideal candidate to 
complement the topologically-protected braiding operations. Moreover, hybrid structures 
involving semiconductor nanostructures and microwave cavities have already been realized 
experimentally [15]. Coupling MBS of a Kitaev chain to a microwave strip-line resonator has 
recently been studied theoretically by Trif and Tserkovnyak [16]. 

We propose to use microwaves to control the coupling between two MBS, which is 
potentially relevant in the context of quantum computation. We will show that photon-assisted 
tunneling has a strong impact on the coupling between adjacent MBS which are separated by 
a short topologically trivial region. In Ref. [17], we presented the main idea using a minimal 
model containing only two MBS coupled via a gapped region with quadratic spectrum. We 
used perturbation theory for weak electron-photon coupling to derive an analytic result for 
the system dynamics, and found that a microwave field can induce Rabi oscillation between 
adjacent MBS. In this paper, we extend that analysis by considering a more realistic model 
and treating the coupling exactly. We explain how coupling between MBS and a microwave 
cavity field emerges in a semiconductor nanowire hosting MBS, and map the corresponding 
Hamiltonian onto a Kitaev chain coupled to a cavity. We solve this model numerically and 
obtain a solution valid at arbitrary electron-photon coupling strength. 

The remainder of this paper is organized as follows. In Section 2 we introduce the 
Hamiltonian for an inhomogeneous Kitaev chain consisting of two topologically nontrivial 
regions separated by a topologically trivial, gapped region. In this situation it is well- 
known that MBS form at the interfaces between the topologically different phases [1]. We 
demonstrate that virtual cotunneling processes mediated by the gapped region couple the 
MBS. The amplitude for these processes decreases exponentially in the length of the gapped 
region, and the decay length diverges as the gap in the topologically trivial region closes. 

The Kitaev model is known to be an effective description for a topological nanowire 
[18]. In Section 3 we study the coupling of such a nanowire to the microwave field inside a 
cavity. We find that the cavity field gives rise to a modulation of the hopping matrix element 
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between neighbouring lattice sites. This enables us to control the properties of the chain with 
the microwave field. 

In Section 4 we discuss a Kitaev model coupled to a microwave field. In the first part of 
that section, we assume that the cavity field can be described classically and investigate the 
dynamics of the driven system. In the second part, we solve the full quantum master equation 
for a damped, quantized photon field. In both cases we find that microwaves can exponentially 
enhance Rabi oscillations in the Majorana sector. They can thus supplement topologically- 
protected braiding operations to gain full single-qubit control which is a necessary step toward 
a universal set of quantum gates for Majorana qubits. 

2. Inhomogeneous Kitaev chain 

To describe a one-dimensional wire which can be brought into either a topologically trivial or 
topologically nontrivial phase by tuning the system parameters, we shall use the Hamiltonian 
of the Kitaev chain [1]. This model captures qualitatively many features of more realistic ID 
models [6, 7]. In fact, it has been shown that there exists an approximate mapping between 
these more realistic ID models and the Kitaev chain Hamiltonian, so their low-energy degrees 
of freedom are identical [18]. We shall extend this mapping to incorporate electron-photon 
coupling in Sec. 3. 

We start with a brief review of some essential properties of an inhomogeneous Kitaev 
chain, containing two topologically nontrivial regions, separated by a short topologically 
trivial, gapped region. We shall derive the effective low-energy Hamiltonian describing the 
overlap between the two MBS adjacent to the central region, and show that such a Hamiltonian 
can in principle be used to implement single-qubit rotations. On a very general level, this 
system can be modeled as a Kitaev chain consisting of iV sites with position-dependent 
parameters, 



Here, /i n denote the onsite energies, t n are the hopping amplitudes, and A n are the p-wave 
pairing strengths. The operator c n (c* ) annihilates (creates) a spinless fermion at lattice site 
n. This bilinear Hamiltonian can be diagonalized exactly for arbitrary parameters and will 
be used below for numerical results. To obtain simple analytical expressions, it is convenient 
to select the simplest set of parameters which generates the desired topological phases. We 
therefore assume Eq. (1) to be of the form 



where the iV-site chain is split into three parts: H L and H R describe the left in < m{) and 
right (n > m 2 ) parts of the chain, and Hq the Nc = m 2 — m\ — 1 sites in the central part 
(mi + 1 < n < m 2 — 1). The sections are coupled by H LC , which connects sites mi and 
mi + 1, and H RC , which connects sites m 2 — 1 and m 2 . 




c n +i + A n c n c n+ i + h.c.) . 



(1) 



Hr = H L + He + H R + Hlc + H RC , 



(2) 
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Figure 1. Band structure of the inhomogeneous Kitaev chain along the y axis: the left and 
right chain segments are in the topologically nontrivial phase and host Majorana bound states 
(MBS) at their edges. The central region is gapped and topologically trivial. Stars denote the 
positions of zero-energy MBS. 



The left and right regions are supposed to be in the topologically nontrivial phase, i.e., 
we assume /i„ = and constant A n = t n > in these regions. Then, the Hamiltonians Hl,r 
become diagonal in a basis of nonlocal Dirac fermions, d n = Im c n+i + iRe c n , 

mi-l N-1 

H L =t L 1 ^2 d l d n, H R = t R 22 d ndn- (3) 

n=l n=m,2 

In the central part of the chain, we choose the parameters as \i n = He, t n = t c > 0, and 
A n = 0. Without loss of generality, we assume jic > 0. Since this parameter choice makes 
the central region topologically trivial, we retain the basis of local Dirac fermions c n , 

m 2 — 1 ^ ma— 2 

H c = ~Hc c « c n - y ^ (4 c n+i + h-c) • (4) 

n=m\+l n=m\+l 

To diagonalize He, it is convenient to extend the central region to a large number of sites 
iVoo 3> N c and then impose periodic boundary conditions. We shall discuss the quality 
of this approximation by comparing it to numerical results below. Hence, we introduce the 
operators c n , where c n = c n+ N x - They are defined to coincide with the original operators c n 
in the central region: — c n for mi + 1 < n < ra 2 — 1. Then, we can diagonalize H c in 
momentum space, 

H c = J2e(k)d{c k , (5) 

k 

where e(k) = —He ~ tc cos(aok) . We introduced the lattice spacing arj, and the momentum 
k is in the first Brillouin zone, k E [— 7r/a , ^/a }, and is quantized in units of 2% / (aoN^) . 
The operators in momentum space are defined by Ck = N^ 1 ^ 2 ^2^=i e~ tha ° n c n . 

If the three segments are not coupled (Hlc = Hrc = 0), our choice of parameters entails 
that the left (right) segment is gapped due to the superconducting pairing and the width of the 
gap is 2th (2tR). The uncoupled Hamiltonian has a four-fold degenerate ground state as there 
are two MBS, 71 and 7 mi , at the ends of the left chain and two MBS, 7 m2 and 77V, at the ends 
of the right chain. In terms of the electron operators, these are given by, 

71 = 21m ci = -z(ci - c\), 
7wi 2Re c mi c mi -\- c m , 
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Figure 2. Splitting of the MBS as a function of chemical potential \ic for different numbers 
of lattice sites Nc in the central region. The crosses indicate the numerically exact solution of 
the inhomogeneous Kitaev chain, see Eq. (1). The solid line represents the splitting calculated 
numerically using periodic boundary conditions. The dashed line shows the approximate 
analytic expression (9) which is valid in the limit /l«c tc- 



1m2 21m c m2 i(c m2 c m2 ), 

7at = 2Re c N = c N + c ] N . (6) 

Now, we introduce tunneling between the side chains and the central chain. Because there are 
no terms beyond nearest-neighbor terms in Eq. (1), only the Majorana states 7 mi and 7 m2 are 
coupled to the central chain [19], 

IT t LC ,~ ~f \ 

n LC — ^7milC mi+ i - C mi+1 J, 

Hrc = -^(c m2 _i + 4 2 -i)7m 2 - (V) 

Note that for our choice of parameters the Dirac fermions d n of the left and right chains are 
not coupled to the remaining system. Therefore, we can discard the operators H LR and use 

Hq = H c + Hlc + H R c. (8) 

The Hamiltonian (8) can easily be solved analytically. We are particularly interested in the 
coupling between j mi and 7 m2 , mediated by the central chain. For fx c < tc, tne central 
chain is gapless and real electrons and holes can tunnel to and from the MBS. In this 
regime, the overlap of the MBS turns out to be very sensitive to system parameters, and the 
boundary conditions of the central region become important. The MBS self-energy acquires 
an imaginary part, indicating a level broadening and thus a finite lifetime of the MBS. In the 
continuum limit, the effect of the metallic central region thus resembles that of a fermionic 
bath which has been investigated in detail in Ref . [20] . 

On the other hand, for jj c > tc, me central chain has a gap of width 2(fi c — t c ). 
The corresponding band structure (for uncoupled chains) is shown in Fig. 1 . In this regime, 
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electron cotunneling is the only process leading to a nonzero overlap of the Majorana levels 
j mi and 7 m2 . This overlap causes a finite level splitting of the MBS, which for \iq tc 
becomes 

e -(A r c -l)arccosh(/i C /*c) 

e = t LC t RC . (9) 

tcV^c/tc) 2 - 1 

Importantly, this splitting decays exponentially with increasing length L c = a N c of 
the central region, which is the reason for the topological protection of the MBS. The 
corresponding decay length £ = a /arccosh(/i C ;/t c -) diverges as (nc/tc ~ l) -1 ^ 2 for 
fJ'C tc- The level splitting (9) can also be found from the exact numerical solution of the 
full Hamiltonian (1). A comparison between the exact numerical result and the approximation 
(9) is shown in Fig. 2. 

In Ref. [17], we calculated the level splitting for a minimal model containing two MBS 
coupled via a gapped region with quadratic spectrum e(k) = k 2 /(2M) — fx . To compare 
the energy splitting (9) with that model, we need to take the continuum limit of the Kitaev 
chain by sending the lattice constant a — > 0, while keeping the total length Lc constant. 
Let us assume that ji c > tc- Near the bottom of the band the spectrum of the uncoupled 
central Kitaev chain is quadratic with effective mass M = l/(agfc) an d chemical potential 
A*o — tc — fJ-c- Then, one finds for small |/x 1 <^C t c , 

e -L c /i 

e w a t LC t RC -, — rr (10) 
2M£ 

where £ = (2M|yU |)~ 1 / 2 is the decay length of the MBS into the gapped system. 

If we focus on the gapped regime > tc), we can derive a simple effective low-energy 
theory of the coupled Majorana modes 7 mi and 7 TO2 by integrating out the central region. This 
leads to an effective retarded interaction between 7 mi and 7 m2 . In the low-energy (long-time) 
limit at energy scales small compared to (fx c — t c ) / nc, we can neglect the retardation and 
obtain the effective low-energy Hamiltonian 
is 

H e S = — 7m 1 7m 2 - (H) 

For the purposes of topologically protected quantum computing, a logical qubit should be 
encoded into four MBS [4], which can in turn be combined into two Dirac fermions, e.g., 
il>L = (li + *7mJ/2 and ip R = (^ m2 + i / y N )/2. Since the topological protection relies on 
a conserved fermion parity, the computational basis should contain states with equal parity, 
e.g., the odd-parity states 

ID = A |o> , 

It) =V&|0) (12) 

where |0) denotes the ground state of the system which is annihilated by ipR an d ?Pl- 
Braiding the MBS 71 and 7 mi (7 m2 and 7^) corresponds to the unitary transformations 
Ui t 2 = exp (±^fa z ) , respectively, where a z is a Pauli matrix in the basis {\\) , \\) }. Similarly, 
braiding j mi and 7 m2 corresponds to U 3 = exp (i7ra x /4). Clearly, these operations are 
insufficient to reach arbitrary points on the Bloch sphere [13]. In order to perform arbitrary 
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single-qubit rotations, one needs to supplement them by topologically unprotected single- 
qubit gates. For instance, if the MBS 7 mi and 7 m2 are subject to a coupling Hamiltonian (1 1) 
for a certain time t, the qubit state will be rotated by the unitary transformation 

U e {t) = exp (-y^J- < 13 ) 

It is easy to show that arbitrary points on the Bloch sphere can now be reached by combining 
the operations U e (t) with a single braiding operation, e.g., U\ [17]. Let us stress again that 
the operation U £ (t) is not topologically protected. In the following, we shall demonstrate 
that the coupling Hamiltonian (11) can be engineered by using the interaction of MBS with 
a microwave cavity field. A protocol to braid MBS using T-junctions of semiconductor 
nanowires was presented in Ref. [18]. By combining this proposal with a microwave cavity, 
as depicted in Fig. 4, it is thus possible to realize arbitrary single-qubit operations. 



3. Coupling of Majorana bound states to photons 

In the proposed solid-state devices [5, 6, 7], MBS exist as quasiparticles consisting of an 
equal- weight superposition of a particle and a hole. Therefore, they may interact with photons 
despite the fact that they are on average chargeless. In order to gain further insight into this 
coupling, we shall use the Hamiltonian of Ref. [7], and couple it to an electric field using 
the minimal-coupling substitution p — )■ p — eA, where p is the momentum operator, e is the 
elementary charge, and A is the vector potential. 

The model Hamiltonian for a semiconductor nanowire along the y-axis reads 

Id 2 d 



(14) 



+ B^l{y)^M + [A ^\(y)^l(y) + h.c. 

<7=t,4- 

Here, p,Q denotes the chemical potential, and u is the strength of the Rashba spin-orbit coupling 
pointing in the z direction. B is a perpendicular magnetic field in the x direction. Finally, the 
induced s-wave pairing strength is denoted by A . It has been shown that, assuming B > A , 
this model contains MBS at the edges for \^ \ < \f B 2 — Aq, whereas it is in the topologically 
trivial phase for |/x | > ^ B 2 — A|. 

Mapping Eq. (14) on a Kitaev chain becomes most transparent in the regime B 3> 
A ,Mw 2 . The band structures in the topologically trivial (nontrivial) phases are depicted 
in Fig. 3. In the nontrivial regime, we can choose /i = and the band gap is proportional 
to A u/B. In the trivial regime, on the other hand, we choose — /i > B and the band gap is 
\/i \—B. One obtains a low-energy theory by retaining the lower branches ^_ w (^— ^)/a/2 
of both spectra [18]. In the nontrivial region, the result is a spinless p-wave superconductor. 
On the other hand, the trivial region becomes a spin-polarized electron system. Both systems 
may be described by Kitaev chains in the respective topological phases. For a fixed lattice 
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Figure 3. Single-particle spectrum corresponding to Eq. (14) for A = 0. Blue lines are 
for the topologically nontrivial phase (/io = 0, B = 2Mu 2 ), green lines are for the trivial 
phase (/io = -2.3M« 2 , B = 2Mu 2 ). For the mapping onto a spinless Kitaev chain, the 
higher-energy parts of the spectrum (dashed lines) are neglected. 



spacing a , the parameters of the homogeneous Kitaev chain 

N N-l 

H K = - ^/ic+c n - - ^2(tc j n c n+1 + Ac n c n+1 + h.c.) (15) 

n=l n=l 

are related to those of the model Hamiltonian H nw as follows, 

1 1 AqU 

» = ^ + B -M^> t= M^ A= B^ - (16) 
For our proposed coupling scheme, the band gap in the nontrivial region should ideally be 
larger than in the trivial regions, so we will assume that |//n|, B 3> A , whereas \no\ — B 
Ao. 

The coupling to a vector potential A(r) in a microwave cavity can now be investigated 
using the minimal coupling substitution —id/df — > —id/df— eA(r) in the Hamiltonian 
(14). We assume that the cavity electric field E = —dA/dt is oriented along the wire 
(y) axis, i.e., the wire axis is perpendicular to the cavity axis, see Fig. 4. Moreover, 
we assume that E{f) is spatially constant along the length of the wire, which is a good 
approximation for current experimental setups [8] and typical microwave wavelengths. In 
that case, A y = (E ims /Cl)(a + a)), where E vms is the root mean square of the cavity electric 
field, a {a)) is the annihilation (creation) operator for the cavity mode, and Vt is its frequency. 
Since A y is position-independent, it commutes with the momentum operator, so the electron- 
photon coupling Hamiltonian becomes 

ieA f 

tfd-ph = ~jf / dy^_(y)d y ^(y). (17) 



Microwave-controlled coupling of Majorana bound states 

A w 20mm 



|7i 



9 




y 



o 



J7m 2 
l7iV 



Figure 4. The microwave cavity is realized as a superconducting stripline. The nanowire 
hosting MBS is located close to one of the maxima of the electric field E(r) and is parallel to 
E. A T-junction in the nanowire allows braiding of two MBS, whereas the microwave field 
allows Rabi oscillations between the two other MBS. 



We can again discretize the spatial integral. For a lattice spacing a , the result reads 

#el-ph = ^2 ( iC n C n+l + h ' C -) °C K a + flt ) C l C n+l + h.C. (18) 

® n n 

In conclusion, an electric field gives rise to a change in the hopping matrix element between 
neighboring lattice sites. The expression (18) holds in both the topological trivial and 
nontrivial regimes. We expect the same type of coupling also at interfaces between different 
phases, and the coupling constant will then depend on the overlap between the wavefunctions 
in both regions and on their spin structure. 

4. Inhomogeneous Kitaev chain in a microwave cavity 

Let us assume that the inhomogeneous Kitaev chain is brought into a driven microwave cavity 
[16]. We consider one cavity mode with frequency and assume that the cavity is driven with 
a frequency VL L , 

F cav (r) = QcJa + ape inLT + aVe"* 1 ^ + # cav , d , (19) 

where r is the time, a is the bosonic operator of the cavity mode, cp represents amplitude and 
phase of the drive, and H cav d contains damping terms, which produce a nonzero line width k. 

We assume that the wavelength of the field A = c/Vt is much longer than the Kitaev 
chain, so we can treat the field as constant along the Kitaev chain. The cavity field acts on the 
electrons forming the Kitaev chain. Thus, we use the Hamiltonian (1) and supplement it by 
electron-photon coupling terms (18). 

For the numerical solution of the system, we represent the Kitaev Hamiltonian (1) as 
H K = \A^U K A, where = ( c\, . . . , cjy, ci, . . . , cm) and T-Lk is a complex 2N x 2N 
matrix. In the absence of coupling between the three regions, the matrix Hk has a fourfold 
degenerate eigenvalue at zero energy reflecting the four MBS (6). The corresponding basis of 
the zero-energy eigenspace is, 

lV>i) = -j= (|1> - |1 + AO) , |W> = ^= (K> + \mi + AO) 

l^iv) = ^ (|A0 + |2A0) , = (|m 2 ) - |m 2 + N)) (20) 



Microwave-controlled coupling of Majorana bound states 



10 



where \n) is the vector with components \n)- = S n j for n E {1, . . . , 2iV}. For the numerical 
simulation, we shall assume that the system is initially prepared in the state \ip mi ). The time 
evolution then leads to Rabi oscillations between the states \ip mi ) and \ip m2 ). 

To apply the proposed coupling mechanism for qubit rotations, the gap in the 
topologically nontrivial left and right regions should exceed the gap in the topologically trivial 
central region, as shown in Fig. 1. Moreover, the photon frequency f2 should be slightly below 
the gap in the central region, and the left and right regions much longer than the central region. 
This choice of parameters ensures that the two MBS in the left region (71 and 7 mi ) and the 
right region (7 m2 and 7^) remain unaffected by the photon field, whereas 7 mi and 7 m2 will 
be coupled. In our numerical simulations, we choose = ji R = 0. In this case, the MBS 
within the left and right regions are mutually uncoupled, so it is sufficient to consider small 
lengths N L)R . 



4.1. Classical microwave field 

Let us now first consider the case where the microwave field can be treated classically. Moving 
into a rotating frame with the drive frequency Q L , the classical approximation corresponds to 
replacing the quantum mechanical operators a, a) by a — > ^yn ph e~ inLT . The Hamiltonian is 
then independent of the cavity resonance frequency f2 and reads 

N 

H = H K - 2/3 v /^cos(fi L r) ^ ( + h.c), (21) 

n=l 

where f3 is the effective electron-photon coupling amplitude and n ph is the number of 
photons in the cavity. Performing a rotating-wave approximation (RWA), we obtain the time- 
independent Hamiltonian 

m,2 — 1 , rri2 — 2 

H RWA = ~ (HC ~ Sl £ ) £ 4 Cn - fJ2 i^n+l + h-C.) 

n=mi + l n=mi+l 



(22) 



2 

Within this approximation, the effective chemical potential of the central region is shifted to a 
new value fi c — > He — i- e - 5 the effective gap of the central region can be tuned by changing 
the microwave frequency Moreover, the coupling strength of MBS to the central region 
can be controlled by n ph which can be changed, for example, by varying the drive strength. 
The Hamiltonian (22) has the same structure as the Hamiltonian (8). The resulting splitting 
between the MBS 7 mi and 7 m2 adjacent to the central region is therefore given by Eq. (9) 
with a shifted chemical potential /jl c — > He — and coupling strength t LC t RC — > /3 2 n ph . 
The low-energy dynamics of the MBS is thus governed by a tunable version of the effective 
Hamiltonian (11), and it is possible to perform the rotations (13) using microwave pulses. 

We will now discuss the validity of the RWA. For this purpose, we initially prepare the 
system in the state \ip mi ). We then start driving at time r = 0, and observe the resulting Rabi 
oscillations of the MBS. In Fig. 5, we plot Rabi oscillations of the MBS 7 mi as a function of 
time, and compare the full numerical solution of the classically driven Hamiltonian (21) to the 
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Figure 5. Comparison between rotating-wave approximation [green line, Eq. (22)] and the 
numerically exact classically driven time evolution [blue line, Eq. (21)]. Here, r is the time, 
and tc the hopping amplitude in the central region. We choose parameters fl^ = Ate and /3 = 
O.ltc, N L = N R = 2, N c = 4, /in = = 0, He = 5t c , A L = A R = t L = t R = 10t c , 
A c = 0. 

corresponding result of Eq. (22). We find that for tc <C < He, both results are in excellent 
agreement. Moreover, the oscillation frequency qualitatively agrees with the approximation 
(9). The deviations from the exact numerical result are of the same order as the deviations 
shown in Fig. 2 and are mostly due to the fact that the analytical result (9) assumed fj, c 3> tc- 
In Fig. 6, we plot the Rabi frequency Qr, determined from the numerically exact solution 
of Eq. (21), as a function of the drive frequency Ql- As expected, we find an exponential 
increase of Qr as VL L approaches the lowest eigenenergies of the central region. In that regime, 
the exact results agree with the predictions from the RWA. If VLl is close to an eigenenergy, VLr 
shows resonant behaviour. In this case, the RWA breaks down. In particular at resonance, i.e., 
if VLl coincides with an eigenenergy of the uncoupled central region, we no longer observe 
clear Rabi oscillations. For stronger electron-photon coupling, resonances also appear if the 
drive frequency matches the halves of the eigenenergies (not shown in the figure), signaling 
two-photon absorption processes which correspond to counter-rotating terms and which are 
neglected in the RWA. 

4.2. Liouville equation in the quantum regime 

The inhomogeneous Kitaev chain coupled to a damped photon mode can be solved 
numerically if we truncate the photon Hilbert space to a finite maximum photon number q c . 
In that case, we can express the system Hamiltonian in the basis \q,n), where q 6 {0, . . . , q c } 
denotes the number of photons, and the single-fermion states \n) for n e {1, . . . , 2A^} were 
defined below Eq. (20). The Hamiltonian matrix thus has dimension 2Nq c . To find Rabi 
oscillations numerically, we solve the Liouville equation governing the time evolution of the 
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Figure 6. Rabi frequency flji as a function of drive frequency for a classical microwave 
field. The yellow shaded region shows the position of the conduction band for a continuous 
central region. The vertical lines denote the eigenenergies of the discrete central region used 
for the simulation. We choose parameters j3 — 0-ltc, Nl — Nr = 2, Nc = 4, /zj, = [lr = 0, 
He = 5tc, Al = A fi = t L =t R = lOto, and A c = 0. 



system density matrix p(r), 
d 



dr 



p(r) = £p(r). 



(23) 



The Liouville superoperator C consists of terms describing the Kitaev chain, the photon field, 
the electron-photon coupling, and the damping, respectively, 

Cp = -i[H K , p] - i[H ph , p] - i[H el>ph , p) + C damping p. (24) 

The Hamiltonian Hk for the Kitaev chain is given in Eq. (1). We again separate the Kitaev 
chain into three segments, with boundaries at 1 < mi < m 2 < N, and choose the parameters 
in such a way that the outer segments are in the topologically nontrivial phase, whereas the 
central chain is topologically trivial. Within each of the segments, the parameters p n and A„ 
are constant. 

The coherently driven cavity mode with resonance frequency Vt is described by the 
Hamiltonian H ph , which contains the photon operators a and a), and a damping term of 
Lindblad form, 

H vh (r) = Slat a + a^(r)e ittLT + ciV* (r)e^ LT , 



£ d ampingP = '^{^pa} - O f ap 



pcJa) 



(25) 



Here, (p(r) encodes the (slowly varying) amplitude and phase of the external drive, VLl is the 
drive frequency, and k is the photon damping rate. Last but not least, the electron-photon 
coupling Hamiltonian has the form derived in Eq. (18), 

N-l 



H, 



el— ph 



-f3(a + a)) (i c t, c n+i + h.c.) 



(26) 
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Figure 7. Damped Rabi oscillations of the Majorana bound states 7 m2 (blue line) and j m3 
(green line) for a quantized photon field with nonzero linewidth (n = 10~ 4 ic)- The cavity 
contains one photon in the initial state and is not driven (ip = 0). The solution was determined 
from the Liouville equation (23). The red line shows the undamped RWA result. We choose 
parameters: fi = 4i c , /3 = 0.1t c , N L = N R = 2, N c = 4, fi L = Mil = 0, [i c = 5t c , 
Al = Ai? = t L = t R = m c , and A c = 0. 

which has been used before in Ref. [16]. We solve the Liouville equation by representing 
the density matrix as a vector with (2Nq c ) 2 components, and the Liouville superoperator as 
a (2Nq c ) 2 x (2Nq c ) 2 matrix acting on this vector. The number of matrix entries scales with 
the fourth power of the Hilbert space dimension, but the matrix is sparse and thus remains 
amenable to a numerical solution. 

The Liouville equation is solved in real time using a fourth-order Runge-Kutta solver. 
We start with an initial fermionic state \ij) mi ). We set the cavity drive to zero and instead 
assume that the initial state contains one photon. Therefore, according to Eq. (20), the full 
system state is |-0j) = (|1, mi) + |l,mi + N)) and the initial density matrix corresponds 
to the pure state pi = \ipi) (ipi\. To observe Rabi oscillations of MBS, we plot in Fig. 7 the 
time-dependent fidelities of the two fermionic states 

Fmifl(T) = Tr [p( T ) (V ® |VW >2 ) (VWx,a|)] , (27) 

where I p h denotes the identity operator in the photon system. The solution of the Liouville 
equation reveals damped Rabi oscillations with a damping rate proportional to k. For weak 
damping, the oscillation frequency is slightly reduced compared with the results for classical 
driving in Sec. 4.1. Moreover, the damping leads to an exponential decay of the amplitude as 
a function of time. 

Within our model, the damping is proportional to k and thus only limited by the photon 
linewidth. Other conceivable damping mechanisms include fluctuations of the photon field, 
disorder in the central region, and electron-electron interactions. However, their investigation 
requires a more detailed study of the actual system hosting the MBS, and is beyond the scope 
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of this paper. 
5. Conclusions 

We have studied an inhomogeneous Kitaev chain consisting of two topologically nontrivial 
regions separated by a topologically trivial, gapped region, and embedded the entire chain in 
a microwave cavity. This extension of the Kitaev Hamiltonian has been shown to provide an 
effective model for a semiconductor nanowire hosting Majorana bound states in the presence 
of a cavity field. We have presented numerical solutions for the cases of classical and 
quantized cavity fields valid at arbitrary coupling strengths. If the microwave frequency 
approaches the band gap of the topologically trivial region, the coupling between the MBS 
adjacent to that region is exponentially enhanced. Switching the photon field on and off can 
then be used to implement controlled rotations of a qubit encoded in Majorana bound states. 
The qubit rotations achievable using this coupling, combined with braiding operations, are 
general enough to allow arbitrary single-qubit gates. 

References 

[1] Kitaev A Y 2001 Physics-Uspekhi 44 131 
[2] BeenakkerC W J 2011 ArXiv:1112.1950 

[3] Leijnse M and Flensberg K 2012 Semicond. Sci. Technol. 27 124003 

[4] Alicea J 2012 Rep. Prog. Phys. 75 076501 

[5] Sato M and Fujimoto S 2009 Phys. Rev. B 79 094504 

[6] Lutchyn R M, Sau J D and Das Sarma S 2010 Phys. Rev. Lett. 105 077001 

[7] Oreg Y, Refael G and von Oppen F 2010 Phys. Rev. Lett. 105 177002 

[8] Mourik V, Zuo K, Frolov S M, Plissard S R, Bakkers E P A M and Kouwenhoven L P 2012 Science 336 
1003 

[9] Rokhinson L P, Liu X and Furdyna J K 2012 Nat. Phys. 8 795 
[10] Deng M T, Yu C L, Huang G Y, Larsson M, Caroff P and XuHQ 2012 Nano Letters 12 6414 
[11] Das A, Ronen Y Most Y, Oreg Y, Heiblum M and Shtrikman H 2012 Nature Physics 8 887 
[12] Majorana E 1937 Nuovo Cimento 14 171 

[13] Nayak C, Simon S H, Stern A, Freedman M and Das Sarma S 2008 Rev. Mod. Phys. 80 1083 
[14] Schoelkopf R J and Girvin S M 2008 Nature 451 664 

[15] Frey T, Leek P J, Beck M, Blais A, Ihn T, Ensslin K and Wallraff A 2012 Phys. Rev. Lett. 108 046807 

[16] Trif M and Tserkovnyak Y 2012 Phys. Rev. Lett. 109 257002 

[17] Schmidt T L, Nunnenkamp A and Bruder C 2013 Phys. Rev. Lett. 110 107006 

[18] Alicea J, Oreg Y Refael G, von Oppen F and Fisher M P A 201 1 Nat. Phys. 7 412 

[19] Bolech C J and Dernier E 2007 Phys. Rev. Lett. 98 237002 

[20] Budich J C, Walter S and Trauzettel B 2012 Phys. Rev. B 85 121405 



